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Abstract 

Soft-gluon resummation within the framework of heavy quarkonium hadroproduc- 
tion is considered. A few selected cases are studied in detail. A sizeable increase of 
the cross sections with respect to the next-to-leading order predictions with central 
factorization/renormalization scale choice can generally be observed. Improvements in 
the dependence of the cross sections on the two scales, especially when they are kept 
equal, are also found. 
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1 Introduction 

In the twilight zone between perturbative and non-perturbative QCD, there hves heavy 
quarkonium physics. Indeed, while the large mass {m ^ Aqcd) of the charm and bottom 
quarks allows for a perturbative evaluation of their production cross sections, still the physics 
of their binding into observable charmonium and bottomonium states unavoidably involves 
low-scale, and therefore non-perturbative, phenomena. 

An understanding of a systematic way to disentangle short- and long-distance effects 
was reached by Bodwin, Braaten and Lepage in ^ . By making use of the Non-Relativistic 
QCD (NRQCD) effective field theory, they were able to provide a framework in which 
heavy quarkonium calculations could be carried out to, in principle, arbitrary high order. 
Non-perturbative soft effects are factored into NRQCD matrix elements, and perturbative 
calculations can provide their coefficient functions. Heavy quarkonium (H) production cross 
sections therefore take the form: 

da{H + X) = J2MQQ[n]+X',fi^){0''[n]y^^\ (1.1) 

n 

In this expression QQ[n] represents a heavy quark-antiquark pair in a given colour and 
spin/orbital/total angular momentum state. Notice that the quantum numbers n might 
differ from the ones of the observable quarkonium H, and at the same time QQ might even 
be in a colour-octet state: soft gluons "hidden" into the non-perturbative matrix elements 
{0^[n]) take care of binding the pair and, at the same time, building up the correct quantum 
numbers, /ia represents the NRQCD factorization scale, separating short- and long-distance 
effects. The relative importance of the various contributions in eq. ( |1 . 1| ) can be estimated 
by using NRQCD velocity scaling rules P], which allow the truncation of the series at any 
given order of accuracy. 

The coefficient functions da{QQ[n] +X) can be evaluated in perturbative QCD. Next-to- 
leading order (NLO) calculations for a wide range of processes in hadron-hadron collisions 
have been presented in ||^. Results for the partonic total cross sections will only depend on 
the scaling variable x = 4m'^/s, s being the partonic centre-of-mass energy. 

It is apparent from the fully analytic results listed in 0] that, when x ^ 1, i.e. when 
the partonic threshold is approached, large NLO contributions can develop and hence spoil 
the convergence of the perturbative series. This calls for an all-order resummation of such 
contributions, which is precisely the goal of this paper. We will work using techniques 
developed in f^, |, ^ §], and more recently applied to heavy-quark production processes 



in I, [11. 

In the following sections we shall first review the structure of the NLO results and the 
resummation formalism, and then present some numerical results. 



2 Next-to-leading Order Results and Resummation 



In ref. [Q results for the total production cross sections &{ij -^ QQ[n] + k) are given 
up to NLO. The colliding partons ij can be gluons, light quark and antiquark, or a quark 
and a gluon. Analogously, the outgoing parton k will be either a gluon or a light quark, and 
results have been given for the QQ pair to be either an S'-state (both scalar and vector, both 
colour singlet and octet) or a ^Pj state. 

For ease of reference let us write down here the cross sections for the processes that 
display the large soft-gluon behaviour. 

The total cross section of the process gg -^ Qg, with Q representing the QQ pair in a 
given state n = ^Sq ' , ^Pq ' , ^P2 ' , reads 

^""[gg^Qg] = <N- Q] (<5(i-x) + ^{ A^otlQ] 5(i-x) 
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In this equation Pgg{x) and Pgg{x) are related to the gluon-gluon Altarelli-Parisi splitting 
vertex: 



Pgg{x) = ICa 



Pgg{x) = 2Ca 



X I - X ,^ , 




1 1 XI ± XI 
1 — X X 

X 1 — X ,^ 


x) 



where the plus-distribution is defined byQ 
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^Notice that in ref. y the resuhs are instead written in terms of p-distributions defined by 

dx{d[x)\pt{x)^ / dxd{x) {t[x)-t(\)\ , 
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'^o[99 ~^ Q] represents the Born cross section for the production of the quarkonium H via 
the intermediate state Q, i.e. according to eq. (|1.1|), 

<b^ - Q] = M99 - Q]{0''iQ)) . (2.5) 

The constants Atot[Q] and the functions fgg[Q]{x) depend, as indicated, on the particular 
QQ[n] state produced. They can be taken, together with aolgg -^ Q], from the results 
published in |^. Finally, n and fip are respectively the renormalization and factorization 
scales. It is worth noticing that the NRQCD factorization scale n\ does not explicitly 
appear, at this order, in the cross section for the production of the QQ pair. It is therefore 
not indicated in the NRQCD matrix element either. 



The presence of very large corrections near a; = 1 is clearly visible in eq. (2J.), in the 
form of (log(l — x)/{l — x))_^_ and (1/(1 — x))_^_ distributions. Their form can be even more 
clearly appreciated by going to Mellin moments space. The Mellin transform is defined by 

fN= [ dx x'^-'fix) (2.6) 

Jo 

and, in the large- A^ limit (corresponding to the x — *> 1 one), eq. (|2.1|) becomes 
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where 7^; = 0.5772... is the Euler constant and the superscript t'^l refers to the colour state 
(singlet or octet) of the QQ pair Q. In this form one can easily see the leading double log, 
due to soft and collinear radiation from the initial light partons, and the subleading single 
logs. The 5c8 in the last line indicates that this term is only present when Q is in a colour- 
octet state, since it is due to soft radiation from a coloured final state. It is worth noting 
that, aside for the constant terms in Atot[Q] and the colour of the final state, this structure 
does not depend on what exactly Q is. 



with p = Aw? / S and S being the hadronic centre-of-mass energy. As a consequence, the constants Atot[Q] 
to be inserted in the NLO cross sections given here differ from the ones pubhshed in [Q , missing now all the 
/^-dependent terms (/3 = -y/1 — p) ■ 



Another process displaying large threshold corrections is qq 
section reads 
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and, as before, 0",^ [gg ^ ^S'j ], Atotp5'| ] and fqq^S{ ]{x) can be found in 
in Mellin moments space. 



(2.9) 
This becomes, 
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Comparing this equation to eq. ( |2.7| ) we can see that the structure is identical, with the 
replacement Ca —>■ Cf for radiation coming from initial-state quarks rather than gluons. 
This cross section can also be seen to be free of an explicit fj,\ dependence. 

Resummation techniques for these large threshold logarithms have been developed in [^, |^, 
1^, 1^ and more recently applied to heavy quark production ||^, [10| . Soft gluon resummation for 
processes also involving non-perturbative matrix elements, namely radiative and semileptonic 
B decays within Heavy Quark Effective Theory, have also been considered in |Tl|, [12|. 

Upon inspection, one can see that the soft limits of the amplitudes projected onto some 
specific QQ[n] state are identical, up to this order, to the ones for open heavy quark produc- 
tion. Such soft limits, derived in |^ , read, for the gluon-gluon channel and for colour-singlet 
and -octet production respectively, 
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where a, b, P, k are, respectively, the momenta of the incoming gluons, the outgoing heavy 
QQ pair, and the emitted soft gluon. A similar expression holds for '^S[ production in qq 
collisions, where we find 
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q and q being the momenta of the two incoming light quarks. 

Making use of these soft matrix elements one can easily reproduce the large- A^ limits 
given above. Moreover, one can then achieve next-to-leading log (NLL) resummation for 
heavy quarkonium production processes by means of the following formula: 

^-[,^- ^ qW^] = ^Q_^[,^- ^ qW] U + ^Q^.fQW]"! A,,,,,^+i (a.(/i), /x, /ip) (2.14) 
the CjjiQ''^^] being the constant, A^-independent terms readable from eqs. (|2.7|) and ( p.lO|) : 
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where the coefficient of the second term is either Ca or Cp for ij = gg or qq respectively, 
and the last term is only present when a colour octet state is produced. 

The resummation function Aij^c.N up to NLL accuracy reads 

Aij,c,iv(as(/i),/i,/iF) = exp|lnA^5ffj^(6o«s(/i)lnA^) + gl^l{boas{n)\nN, jj, fiF)\{'2..lQ) 
The functions g^^^ and g^'^^ resum the LL and NLL terms, respectively. Their explicit form 
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where bo, 61 are the first two coefficients of the QCD /5-function^ 
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It is easy to see that expanding these formulas up to order as, and using eq. ( |2.14| ), one 
recovers the fixed-order results given in eqs. ( |2.7| ) and (|2.10|) . 

Finally, observable hadron-level cross sections will be obtained by convoluting eq. (|2.14|) 
with hadronic parton distribution functions and multiplying by the proper NRQCD non- 
perturbative matrix element. We shall therefore have 

aJ'-^[Q] = F,,^+i(/xf) F„^+i(/xf) a^izj -^ Qg] (0^(Q)) (2.21) 

and an improved hadronic cross section, including the full NLO result plus NLL resummation, 
will be obtained as 

^H,NLO+NLL^g] = a^'''^^[Q] - ((T^'^^^[Q]),3 + ^^'^^^[Q] , (2.22) 

where a^'^^°[Q] is the full NLO resuh as calculated in [|] and {a''''''[ij -^ Qg])a3 is the 
order-a^ truncation of the resummed result, here subtracted to avoid double-counting. 



3 Numerical Results 



Phenomenological cross sections are obtained via inverse Mellin transform of eq. (|2.14 ), to 



be performed numerically. Many problems are to be found at this stage, both of conceptual 
and technical nature. 



^It is worth mentioning the factor of 27r difference between this definition of bo and the one employed in 
ref. i. 
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Figure 1: The cross section for production, in -p-p collisions via the gluon- 
gluon channel, of a fictitious ^^q quarkonium state made up of 50 GeV 
quarks, at 200 and 500 GeV centre-of-mass energy. 



For instance, to take care of the presence of the Landau pole we shall adopt the so- 
called Minimal Prescription introduced in : the integration contour for the inverse Mellin 
transform passes to the right of all singularities in the complex A^-plane, but to the left of 
the Landau pole at A^ = Nj, = exp(l/26oas). 

Technical problems are instead related to the difficulty of performing the required inte- 
grals, due to the resummation function A (a;) being strongly oscillating in x-space close to 
X = 1. The are taken care of by a subtraction procedure similar to the one described in 
Appendix B of ref. . 

In the following we shall present some plots showing the effect of the soft-gluon resumma- 
tion described in the previous section. Hadron level results are obtained with the CTEQ3M 
parton distribution function (PDF) set, unless otherwise stated. NRQCD matrix elements 
for QQ level final states (i.e. not real physical quarkonium states) are fixed, for ^-states, 
at the value number of polarization states x number of colours of the QQ pair. If replaced 
with the real values, in units of GeV^ for ^-states and GeV^ for P-states, the cross sections 
would be in nanobarns. 

We shall first study the renormalization/factorization scale dependence of a fictitious 
heavy quarkonium made up of 50 GeV quarks. Charmonium and bottomonium are very 
close to the non-perturbative region, and the large value of the strong coupling at such low 
scales makes the results less readily readable. 

In fig. |l| resummation effects are shown at 200 and 500 GeV centre-of-mass energies, as a 
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Figure 2: The cross section for production, in pp collisions via the qq 
channel, of a ^S[ state with 50 GeV quarks. 



function of the renormalization/factorization scale /i. Production of a ^Sq state in proton- 
proton coUisions via the gluon-gluon channel is considered here. While the leading order 
result displays the usual monotonical dependence, and the NLO one presents a maximum at 
pretty small values of /i, the resummed result can be seen to be markedly less dependent on 
the arbitrary scale. One can also see resummation effects to be larger closer to the threshold, 
at Vd = 200 GeV, than at 500 GeV, as expected. In both cases a sizeable increase with 
respect to the NLO prediction with /i = 2mQ can be observed. 

Results for other states produced via the gluon-gluon channel, i.e. ^Pq and ^P2, both 
in the colour singlet and octet states, and ^Sq, appear qualitatively similar. These plots 
suggest that a small value of the renormalization/factorization scale should be chosen in the 
NLO calculation in order to obtain a more realistic prediction for the cross section. We can 
also speculate that the value /i = mq, rather than 2mQ, should be chosen as the central one 
for the renormalization and factorization scale, despite the log(4m^//i^) terms appearing in 
the NLO calculation. A variation in the range [mQ/2, 2mQ] would then still span a large 
fraction of the NLO band and give a pretty small uncertainty with the NLO+NLL result. 

Qualitatively similar results are found in the qq channel for production of a '^S'| state. 
Figure |^ shows the dependence of the cross section on the renormalization/factorization scale 
at Vo = 200 and 500 GeV. Again, the NLO+NLL result can be seen to be more stable than 
the fixed-order NLO one. 

One can, of course, also try to vary the renormalization scale n and the factorization 
scale /ip independently. According to the observation made above on the "best" scales 
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Figure 3: Cross section as in fig. ||, but with independent variation of tlie 
renormalization (left plot) and factorization (right plot) scales. Top plots 
are for ^5 = 200 GeV, bottom ones for ^fS = 500 GeV. 



central value, we now fix one of the scales at rriQ (rather than 2?71q) and we vary the other 
in the range [2?tiq/10, 5 x 2mQ]. We can see in fig. |]that the remarkable independence seen 
in fig. |l| is unfortunately at least partially lost, especially when going to large factorization 
scales. Independent scale variations in the qq channel give a similar outcome. These results 



are qualitatively similar to those reported in [^ , where soft-gluon resummation in prompt- 
photon hadropro duct ion is studied, indicating that they are not specific to bound states 
production. They are moreover identical to what can be obtained by making rriQ even as 
large as the top quark mass, i.e. 175 GeV. This suggests that, when studying the overall 
dependence on the renormalization and factorization scales, one should always take care to 
vary them independently. 

We now move to studying real quarkonium states, i.e. made up of charm or bottom 
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Figure 4: The cross sections for production of the bottomonium state ^5q 
via the gluon-gluon channel, at 20 and 50 GeV centre-of-mass energy, as 
functions of the renormalization/factorization scales. 



quarks. In these cases, and especially for charmonium, the not-so-large mass of the heavy 
quark, and hence the relatively large value of the running coupling at these scales, will 
produce less clear-cut results than the ones previously displayed. 

Figure ^ is a scaled-down version of fig. |l], meaning that both the mass of the heavy quark 
and the centre-of-mass energy are scaled down by a factor of 10, to yield 771^ = 5 GeV and 
Vd = 20 and 50 GeV. We can see that, while the behaviour with respect to scale variations 
is the same, even after resumming the soft-gluon logs the cross section does however retain 
a larger dependence than in the rriQ = 50 GeV case, as expected. As before, the NLO-I-NLL 
result predicts larger cross sections than the NLO one with a central scale choice. 

Once again, quahtatively similar results are found for ^Pq and ^P2, both in the colour 

t [81 

singlet and octet states, and for Sq . 

Production of a S[ bottomonium state in qq interactions can also be considered. Fig- 
ure ^ shows the dependence of the cross section on the renormalization/factorization scale 
at y/S = 20 and 50 GeV. At 50 GeV centre-of-mass energy the NLO-I-NLL result can be 
seen to be slightly more stable than the fixed-order NLO one, confirming the trend observed 
in gg collisions. Less is instead to be gained closer to the threshold, at y/S = 20 GeV. 

One can finally also try studying charmonium hadropro duct ion. Figure ^ shows the 
renormalization/factorization scale dependence for production of ^S*,^ and '^S{ charm bound 
states in pp collisions at 20 GeV centre-of-mass energy. The parton distribution function 
set CTEQ41Q (i.e. low Q^, valid down to Q = 0.7 GeV) is used, since very small scales are 
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Figure 5: The cross sections for production, in pp collisions via the qq 
channel, of a ^S[ bottomonium state. 
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Figure 6: The cross section for production, in pp collisions, of ^Sq and 
35f ' charm bound states. The PDF set is CTEQ41Q. 

probed. Moreover, the scales are not varied below /i ~ 0.25 x 2mc ~ 0.75 GeV, already at 
the borderline of both perturbative QCD and the validity range of the PDF set. 

One can see from the plots that the dependence on the renormalization and factorization 
scales is once more lessened when soft gluon resummation is included. However, the im- 
provement is less than in the bottomonium case, as expected, the scale involved here being 
extremely close to the non-perturbative region 
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4 Open Questions 

This paper only provides a first step towards a full understanding of soft-gluon effects in 
heavy quarkonium physics, as a number of questions remain to be addressed. 

The NLO cross sections we have studied here do not present the need for an explicit 
subtraction of infrared singularities to be absorbed into the non-perturbative NRQCD matrix 
elements. Such singularities are however expected to appear in higher orders, as they do 



for instance in the two- loop evaluation of J/ip decay into leptons [|T4|. The resummation 
function A will then have to be modified accordingly, so as to reflect the presence of this 
new factorization scale /xa, in the same way it contains the PDF's factorization scale ^p. 
This will avoid double counting of contributions already included somewhere else. It is also 
conceivable that the interplay of this subtraction with the one of Coulomb terms, which are 
also absorbed, already at the one-loop level, into the NRQCD matrix elements, will have to 
be carefully considered. 

In ref. [Q, where the fixed order NLO cross sections were evaluated, the need for an ex- 
plicit subtraction of infrared singularities, and hence the appearance of the NRQCD factor- 
ization scale /iA, is however in one case already met at order a^. This happens when studying 
the production of a P-wave state via the qq channel. When the emitted gluon becomes soft, 
the ensuing singularity can only be cancelled by adding the corresponding qq — > S-state 
process, in full accordance with the general NRQCD factorization formula ( |1 . 1|) (see also 
section 6 of ref. @] for more details). In this case we find indeed that the soft limit of the 
qq -^ P-state +g process amphtude does not factorize onto qq — > P-state and no gluons, this 
cross section being actually zero, but rather onto the matrix element squared for qq -^ S- 
state. This implies a different framework for achieving soft-gluon resummation, which will 
conceivably be closely interlaced with the structure of the NRQCD matrix elements. 

5 Conclusions 

In this paper we have studied the effect of resumming soft gluons in some selected heavy 
quarkonium hadroproduction processes. The inclusion of resummation effects leads in many 
cases to sizeably larger cross sections than those given by the fixed order NLO result with 
central renormalization/factorization scales. Moreover, improvements in the dependence of 
the cross sections on the renormalization and factorization scales can be observed, especially 
when they are kept equal. We have also pointed out that less marked improvements appear 
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when we vary them independently, leading to an uncertainty band a little larger than what 
would be otherwise obtained. 

A full systematic analysis of soft-gluon resummation in heavy quarkonium processes, with 
phenomenologically relevant results, is in progress, as is the one of the open questions briefly 
described in the previous section. 
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